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Abstract. We construct new integrable coupled systems of = 1 supersymmetric equa- 
tions and present integrable fermionic extensions of the Burgers and Boussinesq equations. 
Existence of infinitely many higher symmetries is demonstrated by the presence of recursion 
operators. Various algebraic methods are applied to the analysis of symmetries, conservation 
laws, recursion operators, and Hamiltonian structures. A fermionic extension of the Burgers 
equation is related with the Burgers flows on associative algebras. A Gardner's deformation 
is found for the bosonic super-field dispersionless Boussinesq equation, and unusual proper- 
ties of a recursion operator for its Hamiltonian symmetries are described. Also, we construct 
a three-parametric supersymmetric system that incorporates the Boussinesq equation with 
dispersion and dissipation but never retracts to it for any values of the parameters. 
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1 Introduction 



In this paper we construct new integrable coupled boson-|-fermion systems of = 1 supersym- 
metric equations and present fermionic extensions for the Burgers and Boussinesq equations. 
The integrability of new systems is established by finding weakly non-local [H 1141 122j recursion 
operators for their symmetry algebras or by describing Gardner's deformations |18| l6l [T2]. We 
find three integrable = 1 supersymmetric analogues of the KdV equation. Also, we relate 
a fermionic extension of the Burgers equation with the Burgers equations on associative algeb- 
ras [lOj. We apply algebraic methods |H] to the study of geometry of supersymmetric PDE, 
and we use the SsToOLS package jSj for the computer algebra system Reduce in practical 
computations. 

First we deal with generalizations of the Burgers equation, which describes the dissipative 
nonlinear evolution of rarified gas. We consider the bosonic super-field version of the Burgers 
equation and construct two infinite sequences of its symmetries; the sequences correspond to 
even and odd 'times' along the flows. Next, we analyze a family of fermionic extensions for 
the Burgers equation itself. We show that if the coupling is zero, then the system at hand 
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reduces to the Burgers super-equation w.r.t. a new field that unites the initial fermionic and 
bosonic components. Otherwise, a new Grassmann independent variable is introduced and the 
fermionic extensions are transformed to the Burgers equation on an associative algebra (see 
also HO]). Also, we observe that an = 2 generalization of the Burgers equation appears as 
a symmetry flow for the Laberge-Mathieu's N = 2 supersymmetric SKdV4 equation The 
diagonal reduction = 0^ of this Burgers system is the A^ = 1 super-field Burgers equation. 

Further, we consider the systems related with the Boussinesq equation, which describes the 
propagation of waves in a weakly nonlinear and weakly dispersive liquid. We construct a Gard- 
ner's deformation |181ll2j of the dispersionless bosonic super-field Boussinesq equation. Thus we 
recursively define its Hamiltonians, which are further transmitted to the boson+fermion repre- 
sentation of the system at hand. Independently, we construct an unusual weakly non-local pil22j 
recursion operator for this system: its differential order is 1, although it proliferates the symme- 
tries of constant order 2. These symmetries are Hamiltonian w.r.t. the previously found chains 
of functionals, and the 'times' along the flows are even and odd variables, respectively. Finally, 
we extend a super-field representation for the 'full' Boussinesq equation with dispersion and 
dissipation to a family of coupled boson+fermion evolutionary super-systems that contain the 
Boussinesq equation but do not retract to it at any values of the parameters. 

The paper is organized as follows. In the introductory part below, we describe three analogues 
of the A^ = 1 supersymmetric KdV equation |15| I16j and construct weakly non-local recursion 
operators for their symmetry algebras, see Example ^ on p. El and Example El in Section [3 
Also, we discuss the properties of a non-local Gardner's deformation ^1 for the A^ = 1 
supersymmetric KdV equation itself. In Section [21 we recall two schemes for generating infinite 
sequences of higher symmetries of the evolutionary super-systems which are contained in the 
experimental database |2I1 and which are the objects studied in this paper. In Section (31 we 
investigate super-field representations and fermionic extensions of the Burgers equation. Then 
in Section El we discuss a Gardner's deformation of the dispersionless Boussinesq equation, and 
we construct a parametric family of super-systems that incorporate the Boussinesq equation 
with dispersion and dissipation. 

1.1 The classification problem 

The motivating idea of this research is the problem of a complete description of A^ = 1 su- 
persymmetric nonlinear scaling- invariant evolutionary equations {ft = (j)^ , bt = 4>''} that admit 
infinitely many local symmetries proliferated by recursion operators; here b{x, t, 0) is the bosonic 
super-field and f{x,t,6) is the fermionic super-field. We denote by 6 the super- variable and we 
put V = Dq + OD^ such that = and [V,!)] = 2D^; here Dq and D X are the derivatives 
w.r.t. 9 and x, respectively (fortunately, the derivative Dg is met very seldom in the text, hence 
no confusion with the operator D occurs). The following axioms suggested by V.V. Sokolov and 
A.S. Sorin were postulated: 



2. All equations are translation invariant and do not depend on the time t explicitly. 

3. None of the evolution equations involves only one field and hence none of the r.h.s. vanishes. 

4. At least one of the right-hand sides in either the evolution equation or its symmetry is 
nonlinear. 

5. At least one equation in a system or at least one component of its symmetry contains 
a fermionic field or the super-derivative V. 

6. The time t and the parameters s along the integral trajectories of the symmetry fields are 
even variables (that is, the parities of b and bt or bs coincide, as well as the parities of /, ft, 
and fs). 



1 



Each equation admits at least one higher symmetry 
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7. The equations are scahng invariant: their right-hand sides are differential polynomials 
homogeneous w.r.t. a set of (half-)integer weights [6] = — |, [x] = —1, [t] < 0, [/], [b] > 0; 
we also assume that the negative weight [s] is (half-)integer. 

Axiom 8 for N > 2 supersymmetric equations that satisfy the above axioms is further de- 
scribed on p. 1141 If a particular equation under study admits symmetry flows with an odd 'time', 
then we use the notation s and the parities of /j, bs are opposite to the parities of / and b, 
respectively. 

Remark 1. From Axiom [S] it follows that the admissible systems are either fermionic exten- 
sions of bosonic systems, or they are > 1 supersymmetric super-field equations and the 
derivations T> are present explicitly. For instance, the superKdV equation @ w.r.t. a fermionic 
super-field f{x,t,6) is = 1 supersymmetric, see ^3 El- On the other hand, an A^ = two- 
component fermionic extension of the Burgers equation that can not be represented as a scalar 
A^ = 1 equation is obtained in Section [3.21 see equation (|16j) on p. 1101 

The first version of the Reduce package SsToOLS by T. Wolf and W. Neun for super- 
calculus allowed to do symmetry investigations of supersymmetric equations. It was used for 
finding scalar fermionic and bosonic super-field and A^ = 1 supersymmetric equations, and for 
description of coupled fermion+fermion, fermion+boson, and boson+boson evolutionary systems 
that satisfy the above axioms; a number of A^ = 2 scalar evolution PDE were also found. The 
bounds < [/], [b] < 5 and > [i] > [s] > —5 were used. The experimental database fI7\ 
contains 1830 equations (the duplication of PDE that appeared owing to possible non-uniqueness 
of the weights is now eliminated) and their 4153 symmetries (plus the translations along x and 
t, and plus the scalings whose number is in fact infinite). 

Remark 2. Axiom [3 together with Axiom El are very restrictive. Indeed, a class of integrable 
systems that do not satisfy these two assumptions is provided by the Gardner's deformations, 
see JHl and (El El- The extended A^ = 1 superKdV^ equation ^ and equation (|27|) on p. El 
give two examples; many other completely integrable extended systems are found in loc. cit. We 
emphasize that no scaling invariance can be recognized for the Gardner's extended systems if 
non-zero values of the deformation parameters £ are fixed. The scaling invariance is restored 
if the weights of the parameters are assumed non-zero: one has [e] = —1 for equation Q and 
[e] = —3 for system (|27|) . Hence we conclude that a classification (see |24j and references therein) 
of the 'symmetry integrable' homogeneous evolution equations may be incomplete, providing 
only zero order terms of the deformations in e. 

Also, we note that the symmetry integrability approach, which was used to fill in the 
database |27j . has revealed a number of systems whose integrability in any sense remains an 
open problem. For example, supersymmetric equation (|33j) admits four symmetries, but this 
knowledge can hardly contribute to constructing a solution of the system at hand. 

Equations ^ and (|33() demonstrate that the database with supersymmetric and coupled 
boson+fermion systems is not exhaustive and, simultaneously, it may contain equations whose 
complete or Lax integrability is uncertain. 

1.2 The N = 1 superKdV equation and its extensions 

The classical integrable supersymmetric evolutionary systems as well as their generalizations 
and reductions are present in the database. 

Example 1 {N = 1 superKdV [151 116|^. Let A^ = 1 and let the weight of the bosonic 
super-field b{x,t,6) be [b] = 1; further, let [t] = —3. Now we scan the cell which is assigned to 
these weights and which is filled in by the runs of SsToOLS. Then we get the potential A" = 1 



4 



A.V. Kiselev and Th. Wolf 



superKdV equation^ 

bt = b^^^ + 3V{b^Vb). (1) 

Indeed, put f{x,t;9) = Vb; then / is the fermionic super-field of weight [/] = | satisfying the 
superKdV equation (jT^EI, see also ^ and [01 Ch.6]) 

ft = fxxx + 3{fVf),. (2) 

The potential super-field equation admits two conserved densities pi = b, p2 = ^b^, which 
are nonlocal w.r.t. equation @; also, potential equation inherits infinitely many fermionic 
conservation laws from superKdV equation . The local conserved densities for are obtained 
by inverting the Miura transformation (see ^l^j) / = X + ^Xx — e^X^X from the Gardner's 
extended = 1 superKdVg equation ^1 EI 



Xt = Xxxx + Hx'Dx),-h''^[{'DxY\ -fe' [x(Px)'J^. (3) 

Recently in ^ it has been observed that the bosonic non-local conserved densities for and 
are obtained by using an appropriate change of coordinates in the potential superKdV^ equation. 
Namely, let us introduce the bosonic super-field (j){x, t, 9; e) such that x = Note that the 
evolution equation for (j) is not in the form of a conserved current, 

(Pt = 4>xxx + 3P(P</> • (t>x) - ^eVx - P(t> ■ €] ■ (4) 

Nevertheless, put tp{x,t,9;e) = exp{e(j)). Then ip{x,t,9;e) satisfies the equation 



,2 



L>Wxx ] h - 



^t = V 

Potential extended superKdV^ equation was implicitly described in p. We note that equa- 
tion (PT| is singular in the deformation parameter e. This is rather unusual with respect to 
the common practice (see and references therein), although singular transformations did 
appear in the context of e-parametric families of integrable equations. 

The Taylor coefficients ^„ of the standard decomposition ^ = ^ • ^n[/] are the re- 

n=l 

quired ^3 nonlocal conserved densities for ((21). The transformation / = f[ip\ defines the 
recurrence relation upon them, and the initial condition is V'l = b. Hence we obtain pi, p2, 
ps = ^b^ + 'D'~^{f'Df), etc. The new densities pn, n > 3 are still non-local w.r.t. potential 
superKdV equation 



1.3 Analogues of the N = 1 super KdV equation 

Now we describe three analogues of superKdV equation ((JJ. These fermionic super-field equa- 
tions are homogeneous w.r.t. the same weights as and they admit infinitely many higher 
symmetries proliferated by weakly non-local O |22l recursion operators. To this end, let the 
weights [/] = I and [t] = —3 be fixed. Then we obtain four evolutionary super symmetric equa- 
tions (namely, ((2) and ((5K,b,c)) that admit higher symmetries fs = (p under the assumption 
[s] > —5. The superKdV equation, see ((21), is the first in this list. We also get the equation 

ft = fxxx + fxT>f. (5a) 



^Throughout this text, the operator "D acts on the succeeding super-field unless stated otherwise explicitly. 
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The recursion operator for equation ()5a() is constructed in Example |21 on p. [71 Third, we obtain 
the two-parametric dispersionless analogue of equation Q: 

ft = a fVU + /? fxVf, a, 13 = const . (5b) 

A computation by Yu. Naumov (Ivanovo State Power University) with SsToOLS demonstrates 
that equation H5b|) admits the weakly non-local recursion operator 

R = afVf ■ V + afU - af^f. ■ - Pf^Vf ■ + • o {fV + Vf) (6) 

and two infinite sequences of symmetries that start from the translations fx and ft- 
The fourth equation for the set of weights [/] = |, [i] = — 3 is 

ft = v{u). m 

It admits the recursion R = fT> — f^T)"^, and also it has an infinite sequence of symmetries 
that starts from the odd weight [s] = —8^. 

Remark 3. The symbols of the evolutionary supersymmetric equations that possess infinitely 
many symmetries are not necessarily constant. For example, equation (|5b|) can not be trans- 
formed to an equation gt = Qxxx + • • • by a differential substitution / = f[g\. The proof is by 
reductio ad absurdum. 



1.4 Remarks 

In this paper, we investigate the geometric properties of the boson-l-fermion systems under the 
additional assumption [/] = [b] (for the primary sets of weights if they are multiply defined). 
From Axioms 01 and 0] on p. [21 it follows that the triangular systems are regarded as trivial and 
therefore their properties are not analyzed. We emphasize that, generally, we can not introduce 
a new anticommuting variable and then unite the two super-fields /, 5 to the fermionic super- 
field (f) = f + Ob such that [c/)] = [f] = [b] — ^ or to the bosonic super-field j3 = b + Of such that 
[P] = [b] = [/] — ^ and such that a scalar equation w.r.t. (j) or (3 holds (there is no contradiction 
with the diagonality assumption because the weights may not be uniquely defined). 

We do not expose now the complete list of supersymmetric boson+fermion systems that 
satisfy the Axioms on p. (21 and such that the weights [/] = [6] coincide. In fact, the symmetries 
for a major part of these equations are proliferated by the recurrence relations (see p. (HI); other 
equations that admit true recursions seem less physically important than the three variants of 
the Burgers and the Boussinesq equations we analyze. 

Yet it is worthy to note some remarkable features of the five systems such that the weight 
[t] = — ^ of the time t is half the weight of the spatial variable x (that is, the equations precede 
the translation invariance). It turned out that these five equations exhibit practically the whole 
variety of properties that superPDE of mathematical physics possess. Let us briefly summarize 
these features. 

Three of the five evolutionary systems are given through 

ft = -afb, bt = b^ + Vf, a = 1,2 A. (7) 

The equations differ by the values a = 1, 2, 4 of the coefficient and demonstrate different 
geometrical properties. The geometry of the a = 2 system is quite extensive: this system 
admits a continuous sequence of symmetries for all (half-)integer weights [s] < — ^, a sequence 
of symmetries such that the parities of the dependent variables are opposite to the parities 
of their flows, four local recursions (one is nilpotent), and three local super-recursions. The 
equation for a = 1 admits fewer structures, and the case a = 4 for equation Q is rather poor. 
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Another equation 

ft = Vb + fb, bt = Vf (8) 

admits local symmetries for all (half-) integer weights [s] < — ^. Equation (jSJ requires intro- 
duction of two layers of nonlocalities assigned to (non)local conservation laws. Four nonlocal 
recursion operators with nonlocal coefficients are then constructed for equation Q. The pro- 
perties of systems 0, ((HJ are considered in details in the succeeding paper 0. 

The fifth system we mention is a super-field representation of the Burgers equation, see 
this system is C-integrable by using the Cole-Hopf substitution. We investigate its properties 
in Section ini 



2 Recursion operators and recurrence relations 

In this section we describe two principally different mechanisms for proliferation of symmetries 
of a PDE. 

2.1 Differential recursion operators 

We consider the (nonlocal) differential recursion operators first. The standard approach OIH] 
to recursion operators is regarding them as symmetries of the linearized equations. The essence 
of the method is the following. The 'phantom variables' (the Cartan forms) that satisfy the 
linearized equation are assigned to all the dependent variables in an equation £; one may think 
that the internal structure of the symmetries is discarded and the (nonlocal) phantom variables 
imitate the (resp., nonlocal components of) symmetries for £. In what follows, the capital 
letters F, B, etc. denote the variables associated with the fields /, b, respectively. Then any 
image TZ = R{ip) of a linear operator R that maps symmetries ip = {fs = F, bg = B} of £ to 
symmetries again is linear w.r.t. the right-hand sides F, B. One easily checks that TZ is then 
the r.h.s. of a symmetry flow 

dsn \bJ 

for the linearized equation hm{£). If the initial equation £ is evolutionary, then the phantom 
variables satisfy the well-known relations 

F[t,s] = 0, B[t^s] = 

that hold by virtue (=) of the equations £ and Lm{£). The method is reproduced literally in 
presence of nonlocalities w whose flows W are described by the corresponding components of 
nonlocal symmetries ip = {fg = F, bs = B, Ws = W). The recursion operators are then defined 
by the triples TZ = (Fg^, Bg^, Wgj^) and generate sequences of nonlocal symmetries. See OEl 
for many examples. 

We finally recall that not each symmetry ip can be extended to a nonlocal flow ip if the set {w} 
of nonlocalities is already defined and, analogously, not all the pairs TZ = {Fg^^, ^g^) generate 
a true recursion TZ. The pairs TZ are therefore called shadows [H] of nonlocal recursion operators. 
The shadows are usually sufficient for standard purposes if they describe the operators that map 
the local components of the flows and whose coefficients are also local. Hence in what follows we 
always set Wg^^ = (that is, we do not find the flows Wg^ that commute with the evolution Wt 
determined by the original system £ and differential substitutions for w). Also, we describe 
the Cartan forms TZ rather than the differential operators R, and we use the term 'recursions' 
instead of the rigorous 'shadows of the generating Cartan forms for nonlocal recursion operators.' 
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Remark 4. The recursion operators considered in this paper are weakly non-local, see I14|. 
Hence one can readily prove the locality of symmetry sequences generated by these recursions by 
using a supersymmetric version of the results in • In the sequel, we use another method for 
the proof of locality. Namely, in Section r4.1.1l we construct a Gardner's integrable deformation [HI 
[T^ [TH] and thus we obtain local Hamiltonian functionals, whence we deduce the locality of the 
corresponding symmetry flows. 

We finally note that a similar method of 'phantom variables' is applied for finding Hamil- 
tonian and symplectic structures for PDE, see [1|. The (nonlocal) Hamiltonian structures for 
supersymmetrizations (|32a|) . (|3.S|) are not extensively studied in this paper. The SsToOLS packa- 
ge is applicable for this investigation as is, since the theory is now transformed to standard 
algorithms of symmetry analysis. 

Example 2. Consider analogue (|5a|) of the superKdV equation ([JJ. We introduce the bosonic 
nonlocality v of weight [v] = 1 such that Dv = f and the fermionic nonlocality w such that 
[w] = I and Vw = (P/)^. In this setting, we obtain the recursion 

n[2] = [Vf-F + + f^.V + ^W). 

The above solution generates the sequence fx^ft^''' of symmetries for equation (|5a|). The 
sequence starts with the translation along x and next contains the equation itself. 

Next, we consider equation ()5b|) and introduce two bosonic nonlocal variables v and w such 
that Vv = f and "Dw = fDf. Then we obtain the nonlocal recursion 

n = a- {VU V - fVfVF + fJF) + (3 ■ {fVf V - f^W) . 

The corresponding operator R is present in ([EJ on p. [51 

Finally, equation (|5bj) is obviously a continuity relation. Therefore, we let v be the bosonic 
variable such that Vv = /; hence we obtain the recursion Fgj^ = fDF — fx V. 

Remark 5. The systems that admit several scaling symmetries and hence are homogeneous 
w.r.t. different weights allow to apply the breadth search method for recursions, which is the 
following. Let a recursion of weight [sr] w.r.t. a particular set of weights for the super-fields /, b 
and the time t be known. Now, recalculate its weight [s^] w.r.t. another set and then find all 
recursion operators of weight [s^]. The list of solutions will incorporate the known recursion 
and, possibly, other operators. Generally, their weights will be different from the weight of the 
original recursion w.r.t. the initial set. Hence we repeat the reasonings for each new operator 
and thus select the weights [sn] such that nontrivial recursions exist. This method is a serious 
instrument for the control of calculations and elimination of errors. We used it while testing the 
second version of the SsToOLS package [Sj. 

The second version of SsToOLS allows to reduce the search for nonlocal recursion operators 
to solving large overdetermined systems of nonlinear algebraic equations for the undetermined 
coefficients which are present in the weight-homogeneous ansatz for Fgj^, Bg^. The algebraic 
systems are then solved by using the program Crack [SHI- The nonlocal variables, which are 
assigned to conservation laws if < 1, were also obtained by SsToOLS straightforwardly using 
the weight homogeneity assumptions. 

Remark 6. The weights of the nonlocal variables constructed by using conserved currents for 
PDE are defined by obvious rules. Clearly, if the weight for a bosonic nonlocality is zero, 
then further assumptions about the maximal power of this variable in any ansatz should be 
made. Within this research we observed that the weights of the new super-fields necessary for 
constructing the recursions are never negative. 
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2.2 Recurrence relations 

Consider an evolution equation E = {ut = 0} and let there be a differential function q[f,h] of 
weight zero w.r.t. an admissible set of weights for 8. Suppose further that the flow Ug^ = q"^ -(j) is 
a symmetry of £ for any n E N; in a typical situation all the flows ipn = Ug" commute with each 
other. Then, instead of an infinite sequence (pn we have just one right-hand side Ug = Q{q) ■ 4> 
of fixed differential order and weight [(p]; here Q is an arbitrary analytic function. In this case, 
we say that the sequence of the Taylor monomials ipn is generated by a recurrence relation. 

We note that the multiplication by q can be a zero-order recursion operator R = q for the 
whole symmetry algebra sym^", otherwise the recurrence relation (pn+i = q ■ fn generates the 
symmetries of £ for the fixed 'seed' flows ipo- The systems that admit recurrence relations for 
their symmetries can possess differential recursion operators as well. 

Example 3. Consider the family of supersymmetric systems 

ft = bVb + fVf, bt = afVb, (9) 

here a € M is arbitrary. We see that equation ® is homogeneous w.r.t. the weights [/] = [b] = ^, 
[t] = —1. The multiplication of the r.h.s. in Q by 6 defines a recurrence relation for infinitely 
many symmetries. Indeed, the flows 

ft = bQ{b)-Vb + f Q{b) ■ Vf, bt = af Q{b) ■ Vb 

commute for all Qs and any constant a. Nevertheless, the operator (q ^) is not a recursion for 
equation Q because it does not map an arbitrary symmetry to a symmetry. 
Further, let a = 1. The system 

ft = bVb + fVf, bt = fVb (Pi) 

admits the local zero-order recursions 

1 _(VbbB + VbfF-VffB\ ^2 _fbxbF + f,bB 
VbfB )' b,bB 



Vbb.,bB + VbbJF - Vbf^fB - Vfb,fB 
13] V -Dbb-JB 



An infinite number of local recursion operators for equation (jy0) is obtained by multiplication 
of TZ^ by 6", n G N. The recursions and TZ^ are nilpotent: (tZ}^ = = (iZf. ^ 



If a = —1, then equation @ also admits infinitely many symmetries that do not originate 
from any recurrence relation because their differential orders grow. 

The recurrence relation ipn^i{ipn, q, n) can depend explicitly on the subscript n; then the 
generators of commuting flows contain the free functional parameters Q{q), Q'{q), etc. 

Example 4. The flows 

'fs{Q)\ _ faf,Qib)+lb,fQ'{b) + 6fb^Q'ib) 



bs{Q)J V abMb)+PUQ'{b) J 

commute for arbitrary functions Q{b) and constants a, /3, 7, (5 G M. Indeed, for any Q{b) 
and S{b) we have 

fAQ)\ (US)\\^ 

br{Q)J'WiS)ji ■ 

The flows defined in (|1U() are also translation and scaling invariant. 
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In the sequel, we investigate the systems that are located on the diagonal [/] = [b] and admit 
infinite sequences of (commuting) symmetries generated by recursion operators; we also analyze 
generalizations of these equations and properties of the new systems. The supersymmetric 
equations whose symmetries are proliferated by the recurrence relations are not discussed in 
this paper. 



3 The Burgers equation 

In this section we investigate three systems related with the Burgers equation. We consider an 
= 1 super-field representation of the Burgers equation and analyze its symmetry properties, 
we relate a fermionic extension of the Burgers equation with the Burgers equation on associative 
algebras, and we indicate an = 2 scalar super-equation whose = 1 diagonal reduction 
{6^ = 6'^) is the bosonic super-field Burgers equation again. 



3.1 Super-field representation for the Burgers equation 

Consider the system 

ft = Vb, bt = b^ + Vf. (11) 

There is a unique set of weights [/] = [6] = ^, [t] = — ^, [x] = — 1 in this case. Hence we conclude 
that the above system precedes the invariance w.r.t. the translation along x. Equation Hll|) 
admits the continuous sequence (|14|) of higher symmetries fs = 4'^,bs = 4>^ at all (half-)integer 
weights [s] < — ^. Also, there is another continuous sequence ()15|) of symmetries for equation (jlH) 
at all (half-)integer weights [s] < — | of the odd 'time' s. 

System Hll|) is obviously reduced to the bosonic super-field Burgers equation 

bx = btt-2bbt, b = b{x,t,9). (12) 

We emphasize that the role of the independent coordinates x and t is reversed w.r.t. the standard 
interpretation of t as the time and x as the spatial variable. The Cole-Hopf substitution b = 
—u~^ut from the heat equation 

Ux = utt 

provides the solution for the bosonic component of ()11() . 

Further, we introduce the bosonic nonlocality w{x, t, 9) of weight [w] = by specifying its 
derivatives, 

Vw = -/, wt = -b; 

the variable w is a potential for both fields / and b. The nonlocality satisfies the potential 
Burgers equation Wx = Wu + such that the formula w = lnu gives the solution; the relation 
/ = —Vw determines the fermionic component in system ()11() . 

Now we extend the set of dependent variables /, 6, and w by the symmetry generators F, B, 
and W that satisfy the respective linearized relations upon the flows of the initial super-fields. 
In this setting, we obtain the recursion 

(Fx-VfF + fxW\ 

'^^'^ = \Bx-VfB + bxW) ^''^^ 



of weight [sr] = — 1. The operator assigned to 71 is 

Vf + fxV-^ 
bxV-^ Dx-Vf 
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The above recursion R is weakly non-local jHEJj that is, each nonlocahty T>~^ is preceded with 
a (shadow [Uj of a nonlocal) symmetry (pa and is followed by the gradient tpa of a conservation 
law: R = local part + ipa ■ 'D~^ o ipa- From j2 it follows that this property is satisfied by all 
recursion operators which are constructed by using one layer of the nonlocal variables assigned 
to conservation laws. The weak non-locality of recursion operators is essentially used in the 
proof of locality of the symmetry hierarchies they generate, see |22j and Remark |3] on p. El 
Recursion ()13() generates two sequences of higher symmetries for system (|11() : 

/A / Vb^-VfVh-Ub \ (U\ (fxx-2Vff,\ 

bj ^ \ VU - {Vf? - b^Vf + bbj ^ ' [bj ^ - 2VfbJ ^ - ^ ' 

The same recursion ()13p produces two infinite sequences of symmetries with the odd parameters s 
for the Burgers equation: 

Vf\ fVU-{Vff - U\ ( fVb-bVf + b, , ^^^^ 



VbJ \T)b,-VfVb-b,fJ ' \bVb- fVf + u- fb^ 

Remark 7. System (lllj) is not a supersymmetric extension of equation (|12|) : it is a represen- 
tation of the bosonic super-field Burgers equation. The flows in p4)) become purely bosonic in 
the coordinates c = T)f, b. The standard recursion for the Burgers equation, see (|20() on p. 1121 
acts 'across' the two sequences in p4|) and maps {ft,bt) i— > {fx,bx)', again, we note that the 
independent coordinates in H12|) are reversed w.r.t. (|19|) . Surprisingly, the flow that succeeds 
the translation along x in (|14|) reappears in (|17j) . 

However, from the above reasonings we profit two sequences of symmetries H15() . which are 
not reduced to the bosonic (x, t)-independent symmetries |10[ § 8.2] of the Burgers equation. 
We finally recall that the Burgers equation 1)12(1 has infinitely many higher symmetries that 
depend explicitly on the base coordinates x, t but exceed the set of Axioms on p. [21 



3.2 Supersymmetric extension of the Burgers equation 

The fermionic extension of the Burgers equation, 

ft = fxx + {bf)x, h = bxx + bbx + a fxf, a G M, (16) 



is a unique extension of the Burgers equation that admits higher symmetries and which was found 
by using SsToOLS. It must be noted that equation (fT6|) contains the unknowns f{x,t), b{x,t), 
and it seems to have nothing to do with a super symmetry. The one-parametric family ()16)) admits 
the symmetries (^'') at all negative integer weights [s] < —1. We also note that system (|16() is 
the Burgers equation itself if the fermionic field f{x,t) is set to zero. 

In what follows, we distinguish the two cases: a = and a G M \ {0} in equation (|16() . We 
claim that the algebraic properties of corresponding extensions p6|) for the Burgers equation 
are then different. Indeed, system (|16() is triangular if the coupling constant a equals zero, see 
Remark |H1 below. 

Obviously, the weights [b] = 1 are [t] = —2 are fixed in system (|16|) . If a = 0, then the 
weight [/] of the fermionic field is arbitrary since the first equation in ()16)) is linear w.r.t. the 
derivatives of /. If a 7^ 0, then the weights of the field f{x,t) and of the dimensional real 
constant a are related by the formula [/] = 1 — ^[a], see Remark |2l on p. Olfor discussion. In 
particular, we have [/] = 1 if a is a true constant w.r.t. the scalings and hence its weight equals 
zero. We conclude that the fields b and / can not be coupled to the super-iield u = b + '&f using 
a super- variable d unless [/] = | and an additional relation unites the weights of "d and a. In 
the sequel, we consider both cases: [/] = | and [/] 7^ |. 
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The fermionic component in H16|) is hnear w.r.t. the field f{x, t), and hence the superposition 
principle is valid for it. The quantity 



is an integral of motion for H16|) . The fermionic variable w{x,t) is assigned to the conserved 
current Dt{f) = D^ifx + bf): we set Wx = f such that w ■ w = and [w] = 0. 

Remark 8. System ()16|) models an unusual physical phenomenon. Assume that at any point 
X G M there are two types of a physical value, the bosonic field with density b{x,t) and the 
"invisible" fermionic field with density w{x,t), and let the dynamics of these two fields be 
described by system ((TB|) (or equation (fT7|) below). Suppose that the initial numeric values 
w{x,0) and b{x,0) of the fermionic and bosonic densities, respectively, coincide ai t = 0. If 
a = 0, then densities will coincide for all t > and the corresponding integral quantities 
w{x, t) dx and b{x, t) dx will be conserved in time. 
If a ^ 0, then the feed-back is switched on in H16|) . A ripple in the fermionic space is the cause 
for the bosonic integral b{x, t) dx to change. Indeed, this quantity is no longer conserved 

unless w{x,t) = const or, generally, unless the condition Wxx{x,t)wx{x,t) dx = const holds 
for all t > 0. We see that the reaction of the fermionic component on the bosonic field depends 
on the incline Wx = f and curvature Wxx = fx but not on the density w{x,t). 

3.2.1 Trivial coupling in (HH): cx = 

Let us suppose that q = 0. Then from H16|) we obtain the system 

wt = Wxx + bwx, h = bxx + bbx- (17) 

System (|17|) appeared in [23] and [2S1 in a different context: both fields w{x,t) and b{x,t) were 
regarded as bosonic, and then Backlund autotransformations |25j and the linearizing substitu- 
tions [2n] were constructed. In Remarkdwe noted that the two-component bosonic system (|17)) 
originates from a flow in 1)14(1 . The field w{x,t) was recognized as fermionic in [3], where 
the Painleve properties of (fT7|) and related systems were investigated. In what follows, the 
field w{x,t) is a fermionic dependent variable. 

Remark 9. The bosonic component b{x, t) in p7|l is linearized by the Cole-Hopf substitution 



where the function q{x,t) is a solution of the heat equation qt = Qxx- From (|18() it follows that 
the Burgers equation is the factor of the heat equation w.r.t. its scaling symmetry. This scheme 
is of general nature and can be used for constructing new equations from homogeneous systems. 
In this paper we do not investigate the relation between the scaling invariance of superPDE and 
their factorizations w.r.t. the scaling symmetries, and we do not study the physical significance 
of these projectivizations and of the new equations. 

Now we reduce ()1()|) for q = to one scalar super-equation w.r.t. a new field which is 
constructed as follows. Let be the new independent super-variable such that 

• its square vanishes, • i? = 0; 

• the variable ■!? anticommutes with the fermionic field / and its derivatives with respect to x; 

• the weight of is [i?] = — ^. 

Recall that the weight of the field / is not fixed if a = 0. Hence we postulate [/] = §• 




/(x, t) dx = const 



b = 2q~'q, 



(18) 
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Further, we define the bosonic field n = 6 + i?/ of weight 1 = [b] = [&] + [/]. Then we obtain 
the Burgers super-field equation again! Indeed, we have 

Ut = Uxx + UUx. (19) 

However, we note that the super-derivative D^-\--d whose square is does not appear in the 
reasonings. Therefore, the weights are motivated by the background concept of supersymmetry, 
but they are not uniquely defined. One could easily fix [/] : = 1 and "i? : = with the same 
equation (|19|) in the end. 

Equation (|19() has the well-known recursion operator 

R = D, + \u+\uxD~^ ll = U, + \uU + \u:,D-^{U). (20) 

The recursion operators for the left- and right-noncommutative ^-Burgers equations have been 
obtained by M. Giirses and A. Karasu (private communication). 

Let us construct an analogue of recursion (|2(J|) for system (|16() if the coupling is a = 0. First 
we introduce the bosonic nonlocality v{x,t) of weight zero by setting = b. The potentials v 
and w satisfy the system 

wt = Wxx + VxWx, vt = Vxx + ^vl- (21) 

Next, we make a technical assumption that the bosonic variables v and V of zero weight appear 
in the recursion for H16|) at most linearly. Then we find out that in this setting there are two 
recursion operators of weight 1. The first operator, 



n 



1 _ , -^wB, + ^bxW - IbxwV + Fx + \UV + \bF - \wbB - \fhV 



[1] 



2Bx + bB + bxV 



is the direct extension of twice the recursion (|2U|) for equation H19|) . Simultaneously, we obtain 
the shadow recursion with nonlocal coefficients, 



7^2 



wBx + bxW + \bxwV + 2Fx + fxV + bF + \wbB + \fhV + 2fB 
[11 " I 



Its differential order is 1; this recursion is not nilpotent. 



3.2.2 Arbitrary coupling in (fTB|l : a G M 

In this subsection we consider the case of an arbitrary non-zero real constant a in system ()16() . 

In view of the preceding subsection, we introduce the independent variable •& that anticom- 
mutes with the fermionic field / and its derivatives w.r.t. x and such that 

• '{}■'{} = that is, the square of "d is now non-zero; 

• the weights of the new variable the dimensional constant a (see Remark I^J, and the 
fermionic field / are related by the formulas 2[??] = [a], [/] = 1 — [d]. 

This time we set [??] = [a] : = 0, [/] : = 1, and we do not consider the square root of based 
on the variable ??. We emphasize that 'd is not an ordinary complex number whose square could 
easily be zero, negative, or positive (and •& would therefore be zero, imaginary, or negative, 
respectively). 

Again, we set u = b+'df; the field m(x, t; •&) is homogeneous of weight 1. Moreover, it satisfies 
the Burgers equation 



Ut = Uxx + UUx, 
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but now equation 1^1)) is an equation on the associative algebra generated by u and its derivatives 
w.r.t. X. This is because for all i ^ j we have^ [^^ij^^j] ~ a 7^ 0; for example, 

uu^ = hh, + + ^ ■ {bj + 6/,) / 66, - ^^fj + ^ • {b.f + bU) = u^u. 

The geometry of equations on associative algebras has been studied recently, see [^, by using 
standard notions and computational algorithms. 

Obviously, the fermionic nonlocality w such that Wx = f is indifferent w.r.t. the value of the 
coupling constant a. We find out that system (|16() always admits the conservation law that 
potentiates the bosonic variable 6. We thus set 

Vx = b+ \ fwa, vt = bx + \b^ + \ fxwa + \ fbwa. 

The weight of the nonlocality v is zero. Surprisingly, the variable v satisfies the same potential 
Burgers equation as the nonlocality see equation (f2T|l . 

Wt = Wxx + VxWx, Vt = Vxx + \vl. (HH) 

We conclude that equation potentiates system ((TB|) for all a = • although the algebraic 
nature of the Burgers equation ()19|) w.r.t. the field n = 6 + ?9/ is radically different from 
equation (^I)) that describes the Burgers flow on the associative algebra. 

In the nonlocal setting {/, i«} + {6, there are two generalizations of recursion ()2U|). The 
first recursion of weight 1 for equation ()11() is 

^1 ^( bxW + 2Fx + fxV + yxwWa + bF + fB \ 
[^1 ~ \2Bx±bB_+bxy_ + ^bxwWa + fxWa - fFa) ' 

here we underline the component that corresponds to (\'20\\. The second extension of weight 1 is 

n^ll = -wBx - \bxwV + \fxwWa - \wfbWa - \wfFa - \whB - \fbV - fB, 
7^J^j' = 2Bx + bB + bxV - wFxU + fxWa + IfxwV + IfbWa - \wbFa 
- \wfbVa - IwfBa. 

Remark 10. In the previous reasonings we treated system ()16|) as an = equation that 
involves the fermionic field but does not contain the super-derivative P. Now we enlarge the 
(x, t, f, 6) jet space with the anticommuting independent variable 6 and the derivatives of / 
and 6 w.r.t. 9. We have = Dx, and the unknown functions become the super-fields f{x,t,6) 
and b{x,t,6). Physically speaking, we permit the consideration of conservation laws at half- 
integer weights for H16|) and ij^ff)) . We discover that there are many nonlocal conservation laws 
for equation (|16|) : for example, we obtain the 'square roots' of the variables v and w. We 
conjecture that there are infinitely many = 1 conservation laws for equation (|11() . Also, there 
are many recursions that involve the nonlocalities assigned to the new conservation laws and 
which are nilpotent if a = 0. This situation is analogous to the scheme that generates non-local 
Hamiltonians for the A^ = 1 superKdV equation (jSJ, see p.EJand [T| ITT]. 

3.3 N = 2 supersymmetric Burgers equation 

Let us recall that Dgi and Dx denote the derivatives w.r.t. the independent coordinates 0* and x, 
respectively, while are the super-derivations such that T)f = Dx for any i. 

We now admit that Axiom 8 was used when constructing the evolutionary super-systems: 

^One could use the notation h instead of a. 
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8. Each of the super-derivatives Vi = D^i + O^D^, i = 1, . . ., N , occurs at least once in the 
r.h.s. of the evolutionary system if > 2. 

In the database [^Tj there is a scalar, third order N = 2 supersymmetrization of the Burgers 
equation; it is 

bt = ViV2b^ + bb^, b = b{x,t,9\e^), V, = Dg,+9'D^, i = l,2. (22) 

Equation (|22p is reduced to the second order Burgers super-field equation bt = b^x + bb^ on the 
super-diagonal 6^ = 9"^, here b = b{t, x, 9^, 9^). 

Let us expand the super-field b{x,t;9^ ,9^) in 9^ and 9^: 

b = (3{x, t) + 9^i{x, t) + 9'^T]{x, t) + 9^9'^-f{x, t). 

Hence from equation (|22|) we obtain the system for the components of b: 

Pt = -Ix + Pfix, it = Vxx + {POx, Vt = -ixx + {f3v)x, 

It = f3xxx + {P7)x - {iv)x- (1221) 

We see that the second and third equations in system ((^^ are Burgers-type, that is, they 
contain the dissipative terms and the remaining parts are total derivatives. The fourth equation 
in (|^^ describing the evolution of 7 is of KdV-type: the dispersion and two divergent terms 
are present in it. 

The KdV nature of the N = 2 supersymmetric Burgers equation (|22j) is not occasional. 
Indeed, equation is a symmetry of the N = 2 supersymmetric SKdV4 equation JH] 

bs = -bxxx + \{bViV2b)^ + l{VibV2b)^ + IbHx. (23) 

Reciprocally, the SKdV4 equation (|23jl is a higher symmetry of the Burgers equation (|22|) . and 
their bi-Hamiltonian hierarchy is of the form 

bx ^ equation 1-^ equation i-^^ • • • . 

The two equations share the recursion operator |13| I17j . It must be noted that the relation 
between the Laberge-Mathieu's N = 2 SKdV4 equation and the N = 2 Burgers system was not 
indicated in loc. cit. 



4 The Boussinesq equation 

In this section we describe a super-field representation of the dispersionless Boussinesq equa- 
tion; the system at hand admits two infinite sequences of commuting symmetries of constant 
differential order 2 which are generated by a weakly non-local recursion operator of differential 
order 1. Also, we extend the Boussinesq system with dispersion and dissipation to its three- 
parametric analogue that does not retract to it for any values of the parameters. 



4.1 The dispersionless Boussinesq equation 

We consider the two-component system 

ft = bVb, bt = Vfx. (24) 
System (|24() is a super-representation of the dispersionless Boussinesq equation 

btt = \{b')xx, (25) 

here 6(x, t, 9) is the bosonic super-field. In what follows, we construct a Gardner's deformation 
for the hydrodynamic representation of the dispersionless Boussinesq equation (|25|) . Next, we 
transmit the properties of Hamiltonian symmetries for equation (|25jl onto its supersymmetric 
representation (^l)) . 
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4.1.1 The Gardner's deformation of the dispersionless Boussinesq equation 

The bosonic two-component form of ()25() is 

bt = Cx, ct = hbx- 



(26) 



Here we obviously have c{x,t;6) = X'/(x,t;6'). System is a super-field equation of hydro- 
dynamic type. 

Remark 11. Consider the dispersionless Boussinesq equation (|26|) with b{x,t) and c{x,t). The 
number of independent variables in it coincides with the number of unknown functions and equals 
two. Therefore, the system at hand is linearized [22 by using the hodograph transformation 
b{x,t), c{x,t) I— c), t{b,c). Indeed, we obtain the linear autonomous system 



Xh 

b ' 



ESI) 



here c is the new time and b is the new spatial variable. The solution of the hydrodynamic type 
system is expressed with the Airy function. 



Now we construct a Gardner's deformation pi 1121 IT^ for system ()26|). We have 

n{b,c), (EniO 





D,, 



5/5b 
5 /5c 



{x,t) 



where the density of the Hamiltonian Ti is H = ^b'^ + hc^. We proceed with a deformation 



w 





D, 



ix,t) 



of equation (|26r|) . We assume that the deformations of the dependent variables are 
b = w'^+ eWi{w) + e'^W2{w) -\ , c = w'^ + eQi{w) + e'^Q.2{w) H , 

and we extend the Hamiltonian functional such that its density is 
H = H{w) + eHi{w) + e^H2{w) + ■■■ . 

We also expand the fields w^{x,t;e) and 'uP'{x,t;e) in e: 



(27) 

(28a) 
(28b) 



w\x,t]£) = ^ 



Wi. e 



1,2. 



fc=0 



Recall that equation H26|) is in the divergent form. Therefore, the Taylor coefficients w\. are 
termwise conserved, and from (|28aj) we get the initial conditions Wq = b, Wq = c and the 
recurrence relations for the conserved densities. 

Now we truncate expansions (|28a|) . ()28b|) to polynomials of sufficiently large degrees. By 
using the homogeneity reasonings, in view of [e] = —3, we then reduce the deformation problem 
to an algebraic system for the undetermined coefficients in these expansions. We solve the 
algebraic system by using the program Crack [201 ™d finally obtain the deformation 



b = + ew^up' , 



w 



2\3 
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The initial terms in the two sequences of the conserved densities are 
ifg = b, wl = —be, W2 = 2bc^ + i^'*, . . . , 

wl = c, w\ = -(? - i6^ wl = \(? + |6^c, .... (29) 

Taking into account that the Hamiltonian operator ( ^ maps gradients of conservation laws 
for equation to its symmetries, see Pj, we thus conclude that describes two infinite 
sequences of Hamiltonians for l|26|) . They determine infinitely many contact symmetry flows, 
which are local w.r.t. b and c. Moreover, each flow lies in the image of by construction, 
and hence they induce local symmetry transformations of the super- variable /(x, t; 0) such that 
c = P/. We emphasize that we did not even need a recursion operator for H26|) to obtain the 
contact symmetry flows and prove their locality. 

Remark 12. The problem of constructing integrable deformations for homogeneous (super-) 
PDE is another practical application of the Crack solver for large overdetermined systems 
of nonlinear algebraic equation. This application has not been previously considered within the 
framework of 01131271. 

4.1.2 Supersymmetric representation of equation (|25|) 

Now we discuss the Hamiltonian properties of system (|24() . It is homogeneous w.r.t. multiply 
defined weights; we let the primary set be [/] = [6] = 1, [t] = — 1^, [x] = — 1. With respect to 
this set, system (|24|) admits two infinite sequences of Hamiltonian symmetry flows of unbounded 
weights —1, —1^, —3, . . ., and constant differential order 2. These two sequences start 

with the flows 

/A /Vbb^ + VffA //A /VfVbb+y^b^\ 

j>J ^ [vf, b + vfbj^ ' \bj ^ \vu vf + Ib^by ^ ■ ^ ^ 

The skew-adjoint operator A = (_p is a Hamiltonian structure for the symmetries in (|3fl|) . 
and, by Section [4. 1.11 all the flows possess the Hamiltonian functionals. There are two Casimirs 
H^Q^^ = b, H^^^ = Vf for equation (EH). From ^ we obtain the Hamiltonians with densities 

Hi'^ = bVf, Hi'^ = ^b^ + lbiVf)^ 

iff) = i(p/)2 + i63, H^') = i(^vff + ib'vf, .... m) 

Let us construct the nonlocal recursion that maps the symmetries in ()3U() . To this end, we 
introduce the nonlocality w such that wt = b^/2 and Dw = f. We further let the variable v be 
such that vt = T>f, Vx = b. Then we obtain the nonlocal recursion 

^ ^/VbbV+'^b^VV + lVfF + lfxW\ 
\VfxV + ^bVF+lVfB + lbxW J 



of differential order 1. Indeed, the order of the flow 6^^+! equals the order of /s- plus 1 owing to 
the presence of the underlined differential operator V = Dg + 9Dx in TZ^. 

Proposition. Recursion (|31j) proliferates the symmetries of constant differential order 2 for the 
dispersionless Boussinesq super-field equation ()24() . The initial terms of the symmetry sequences 
are given in (|3(H) . and their Hamiltonians are 

The assertion follows from the results of Section [4. 1.11 

We note that Hamiltonian symmetries (|3U() are an example of infinitely many flows that 
are not obtained by a recurrence multiplication scheme, see p. |H1 for definition, although their 
differential order is constant. Indeed, they are contact in the coordinates b, c and are of second 
order w.r.t. / and b. 
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4.2 The Boussinesq equation with dispersion and dissipation 

Representation ()24() of the dispersionless Boussinesq super-field equation is embedded in the 
one-parametric family of supersymmetric systems 

ft = b'Db + 'Dbxx-afxx, bt = Vf^ + abxx, a G M. (32a) 

By definition, put c = T>f . Then from l|32a() we get the bosonic Boussinesq system (see, 
e.g., 0111111) 

ct = bbx + bxxx - acxx, h = Cx + abxx, a G M. (32b) 

The two systems (|32|) are homogeneous w.r.t. the weights [b] = 2, [/] = 2^, [c] = 3, and [t] = —2. 
We see that equation ()32a() has no nontrivial bosonic limit. Indeed, one can not set / = such 
that the fermionic equation remains consistent unless b = const. 

If Q = 0, then (|32|) is the Boussinesq equation with dispersion. The terms involving a describe 
the dissipation. There is a well-known recursion operator of weight [s/jj = —3 for the Boussinesq 
equation without dissipation The recursion for the fermionic component in H32a|) is then, 
roughly speaking, the component corresponding to c in the recursion for ()32b|) conjugated by D. 

If a is non-zero, then system ()32b|) is not reduced to a scalar fourth order super-field equation. 
In this case, the system is translation invariant and for all a G M it admits symmetries of weights 
[s] = - {3k + I ± I) , k € k > 0. 

4.3 The multi-parametric Boussinesq-type equation 

Finally, we construct the Boussinesq-type system using representation (|32a() for the Boussinesq 
equation. Prom the database j27j we obtain the system 

ft = aPfb - a-fbVb - -/"^Vbxx - P^fxx, 

bt = al3b^ + (3^Vfx + /3^bxx, (33) 

here a, /3, 7 G M. This system is an analogue of Boussinesq equation (|32|) but does not retract 
to it for any value of the parameters. Similarly to equation H32a() . system H33() does not have 
a nontrivial bosonic limit at / = 0. 

The translation invariance of equation (|33j) is obvious; we recall that [x] = —1 and [t] = —2. 
Also, equation H33|) admits the symmetry 

/, = -Vbbxxl^ + VbxbxJ^ + VbxVfPj^ - VfxVbp^^ - Vfxfp^j + Vffx(3^j 

-bxxfh^ + bxfx(3i^, 
bs = -VbfxP^j + VbxfPh + -DbxVbp-f^ + fxfp^ 

such that [s] = —4, and there is the symmetry 

/, = VbfxP^' - VbxfUl^ - -DbxVbj' - bh'' - {Vff[i^-i - fxff3'j - 2Vfbxh\ 
b-s = VbbxfSj^ + VfVbpS + VffP^^ + bxf/3h, 

here the weight of the odd parameter s is [s] = —3^. No more symmetries exist for equation (|33() 
if [s] > —7 and if the parity of s is arbitrary. Quite strangely, there are no conservation laws for 
equation (|33jl within the weights not greater than 11; hence no nonlocalities were constructed 
and no recursion is currently known for system H33() . 
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5 Conclusion 

In this paper, we investigated the integrabihty of fermionic extensions and supersymmetric 
generahzations for the KdV, Burgers, and Boussinesq equations. Recursion operators for their 
symmetry algebras were obtained. Also, we analyzed the properties of Gardner's deformations 
for the dispersionless Boussinesq equation ()26() and = 1 supersymmetric KdV equation 

The boson+fermion super-field representations ()11() . H24|l . and (jSH) for the Burgers and 
Boussinesq equations are remarkable by themselves. Indeed, the roles of the independent vari- 
ables X and t are swapped in system ()11() : x is the time and t is the spatial coordinate in 
equation (|12j) . The dispersionless Boussinesq system (|24|) admits two recursive sequences (|30|) 
of Hamiltonian symmetries whose differential orders is constant, while the differential order of 
nonlocal recursion (|3H) is strictly positive. The Boussinesq-type system (|33j) is multi-parametric 
and contains the Boussinesq equation with dispersion as a component but can not be reduced 
to it at any values of the parameters. We finally note that all the Boussinesq-type systems (|24|1. 
H32a() . ()33() . as well as representation 1)11(1 for the Burgers equation, do not have bosonic limits 
at / = 0. 

Using fermionic extension (|16|) of the Burgers equation, we conclude that the 'direct A^ > 1 
supersymmetrization' ^Hl based on replacing the derivatives with super-derivatives in a PDE is 
not a unique way to obtain its generalizations. Indeed, the new systems can contain the fermionic 
fields but no super-derivatives. Hence we indicate three possible types of these generalizations. 

1. The new super-fields combine the old fields with the new added components and satisfy 
manifestly A^ > 1 supersymmetric equations, see the superKdV equation @, for example. 

2. The evolution of new fermionic fields f{x, t) is coupled with the original equation, see ()16|) . 
Then, new independent (Grassmann, or super-) variables -i? are introduced such that the 
extended system is reduced to a smaller equation. The geometric structures of the exten- 
sion are then inherited from the final system by routine expansions in the new variables t?, 
see (fT9|l . 

3. Similarly to the previous case, the new field is constructed by using a new independent 
variable that anticommutes with the fermionic function but not with itself. Then the 
resulting system is an evolution equation on an associative algebra [10] , see equation (^1)) 
on p.d 

We conclude that the equations on associative algebras present a nontrivial way to generalize 
coupled fermion+boson systems. 
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